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Single-electron ionization and excitation cross sections as well as cross sections for excitation 
into the first excited p state of the alkali metal atoms Li(2s), Na(3s) and K(4s) colliding with 
antiprotons and protons were calculated using a time-dependent channel-coupling approach. For 
antiprotons an impact-energy range from 0.25 to 1000 keV and for protons from 2 to 1000 keV was 
considered. The target atoms are treated as effective one-electron systems using a model potential. 
The results are compared with theoretical and experimental data from literature and calculated 
cross sections for antiproton-hydrogen collisions. For proton collisions a good overall agreement is 
found which confirms the present numerical approach, whereas discrepancies are found between the 
present antiproton cross sections and those calculated by Stary et al., J.Phys. B 23, 263 (1990). 
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I. INTRODUCTION 



Collisions with alkali metal atoms as targets have 
been studied in numerous experimental and theoretical 
works over many years. Among these studies a num- 
ber of efforts deal with proton-alkali metal atom colli- 
sions 0, B II EL S B EL B B E3, El Ei El El EE p 

and a smaller number of attempts address collisions in- 
cluding antiprotons as projectiles [13, EH ■ One reason 
for the attractivity of alkali metal atoms is that they 
are relatively easy to access experimentally as well as 
theoretically which opens up the possibility for detailed 
comparisons. The given shell structure of the alkali met- 
als suggests in a theoretical description the application 
of a quasi-one-electron model for the outermost loosely 
bound electron. The electron is then described by means 
of a model potential formed by the Coulomb potential of 
the nucleus and an effective potential representing the 
frozen inner-shell electrons. In particular, Li and Na 
atoms colliding with protons and electrons have been 
in the focus of the investigations so far. Whereas the 
literature on alkali metal atoms collisions dealing with 
antiprotons as projectile is still sparse compared to the 
treatment of protons and electrons. In order to obtain 
cross sections for ionization of Li by antiproton impact a 
continuum-distorted- wave eikonal-initial-state model has 
been used by McCartney et al. Furthermore, an 

optical-potential description of collisions of antiprotons 
with Li and Na has been provided by Stary et al. [l7l ]. 
No experimental data are available for the considered 
antiproton-alkali metal atom collision systems yet due 
to the lack of appropriate low-energy antiproton sources. 
This may also be the reason for the relatively small in- 
terest in antiproton-alkali metal collisions compared to 
their proton counterparts until now. However, the up- 
coming Facility for Antiproton and Ion Research (FAIR) 
with its incorporated Facility for Low energy Antiproton 
and Ion Research (FLAIR) [19( will provide the necessary 
experimental conditions in the near future and is there- 
fore expected to attract further attention to the field of 
antiproton collisions. 



The primary motivation of this work is to shed more 
light on the antiproton-alkali metal collision systems and 
to provide a consistent data base for Li(2s), Na(3s) and 
K(4s) atom collisions with antiprotons and protons in a 
large energy range. It starts at low energies E = 0.25 keV 
where the collision processes depend considerably on the 
projectile and ranges up to high energies E = 1000 keV 
where the antiproton and proton collision systems are 
supposed to show the same behavior due to the expected 
applicability of the first Born approximation. The calcu- 
lations for collisions with proton projectiles are consid- 
ered to be valuable in two aspects. On the one hand, 
the proton results - especially for Li targets - can be 
compared in detail with literature values. This way the 
proton results can be utilized in order to test the present 
method and its implementation which is the same for 
protons and antiprotons. Furthermore, new theoretical 
ionization and excitation cross sections for proton col- 
lisions - especially for K targets - are provided which 
to the authors' knowledge were not fully known in the 
energy range considered here. 

Besides the obvious similarities of protons and antipro- 
tons as projectiles they differ mostly in their capture be- 
havior. First, only antiprotons can annihilate with pro- 
tons of the atomic nucleus. Since it is known that the pro- 
cess of annihilation is only likely to occur at very low en- 
ergies [20] it is not included in this investigation. Second, 
in the case of proton collisions electron capture by the 
projectile from the target atom is possible. This process 
plays a dominant role for low-energy collisions. Hence, 
a two-center approach appears to be most promising for 
low-energy proton collisions. However, at low energies 
the present calculations concentrate on antiproton colli- 
sions only. Therefore, a basis expansion which is centered 
solely on the target alkali metal atom for both antipro- 
ton and proton projectiles is used. Thereby, limitations 
pertinent to a molecular approach at high energies are 
avoided. Furthermore, the same method can be used for 
antiproton and proton collisions which may be confirmed 
by a detailed comparison of present proton results with 
literature data. A detailed analysis of the electron cap- 
ture process for proton scattering, however, lies beyond 
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the scope of this work. 

The paper is organized as follows: Sec. [TT] explains 
how the alkali metal atoms are described and reports on 
the computational approach. Sec. IIIII considers the con- 
vergence behavior of the present antiproton and proton 
results. Subsequently, the calculated cross sections are 
presented and compared to literature data. Finally, the 
present results for antiproton-alkali metal collisions are 
discussed and a comparison with a hydrogen atom as tar- 
get is made. Sec. IIVI concludes on the present findings. 
Atomic units are used unless it is otherwise stated. 



II. METHOD 

In this work the target atoms are treated as effective 
one-electron atoms. The valence electron is exposed to 
a model potential V mo d suitable for alkali metal atoms 
which describes its interaction with the nucleus as well 
as with the remaining core electrons. Additionally, core 
polarization effects are included. The employed model 
potential was proposed by Klapisch plj . The used poten- 
tial parameters are given in [22| . The effect of the spin- 
orbit coupling is not included in the present approach. 
The energies of the states with principal quantum num- 
ber n < 6 for the alkali metal atoms Li, Na and K which 
were achieved with this approach are given in table U to- 
gether with compiled values of the NIST data bank [23| • 
In the case of energy level splitting due to spin-orbit cou- 
pling the present energies are compared to the lower lying 
reference energies. The largest relative energy splittings 
of the reference data due to the spin-orbit coupling of 
the energies given in tabled are 0.002 %, 0.1 % and 0.4 % 
for the energetically lowest lying p states of Li, Na and 
K, respectively. Particularly for Li there is a very good 
agreement with the data given by NIST. But also for the 
other two atoms the deviation from the literature values 
remains at maximum around one per cent. 

In order to describe the collision process the relative 
motion of the heavy particles is approximated by a clas- 
sically trajectory (CTA) also referred to as the impact- 
parameter representation. The projectile is assumed to 
move on a classical rectilinear trajectory with a constant 
velocity v parallel to the z axis. The internuclear dis- 
tance vector R is given by R = b + vt, where b is the 
impact-parameter vector along the x axis and t the time. 

The time-dependent Schrodinger equation 

i-V(r,t)= (H +V int (r,R(t))) tt(r,t) (1) 

of the target atom interacting with the projectile is 
solved. The atomic Hamiltonian of the target atom is 
defined as 

H = -iv 2 +V mod (2) 
and the time-dependent interaction between the projec- 



n I 


Li calc. 


Na calc. 


K calc. 




Li ref. [23] 


Na ref. [23] 


K ref. [23] 


2s 


-0.198477 
-0.198142 






2p 


-0.130482 
-0.130236 






3s 


-0.074362 
-0.074182 


-0.189163 
-0.188858 




3p 


-0.057364 
-0.057236 


-0.111760 
-0.111600 




3d 


-0.055605 


-0.056071 


-0.061596 




-0.055606 


-0.055937 


-0.061397 


4s 


-0.038694 


-0.071754 


-0.160105 




-0.038615 


-0.071578 


-0.159516 


4p 


-0.032036 


-0.051075 


-0.100434 




-0.031975 


-0.050951 


-0.100352 


4d 


-0.031274 


-0.031531 


-0.034954 




-0.031274 


-0.031442 


-0.034686 


4f 


-0.031254 


-0.031267 


-0.031337 




-0.031243 


-0.031268 


-0.031357 


5s 


-0.023677 


-0.037656 


-0.064121 




-0.023637 


-0.037584 


-0.063713 


5p 


-0.020408 


-0.029265 


-0.047187 




-0.020374 


-0.029202 


-0.046969 


5d 


-0.020013 


-0.020160 


-0.022158 




-0.020013 


-0.020106 


-0.021983 


5f 


-0.020002 


-0.020010 


-0.020048 




-0.019969 


-0.020011 


-0.020062 



Table I: Calculated binding energies (Hartree) for Li, Na and 
K using a Klapisch-model potential. The reference data is 
taken from the NIST data tables [23l ]. In the case of energy 
level splitting due to spin-orbit coupling only the energetically 
lower lying reference energy is given. 

tile with the charge Z p and the target atom as 

where r is the spacial coordinate of the explicitly treated 
valence electron. 

The total wavefunction ^(r, t) is expanded as 

*(r,t) = ^2 c nim(t) <?W(r) exp(-ze„ z i) (4) 

nlm 

using the expansion coefficients c n i m (t). Here, the <f> n i m 
are eigenstates of the Hamiltonian Ho with the energy e n i 
obtained with the model potential. The n, I, m are the 
principal one-electron quantum number, angular momen- 
tum, and its projection on the z axis, respectively. The 
4>nim can be further expanded as 

•Wr) = g n i(r) I [(-l) m YT(Sl) + y ; m (f2)] , 

a/2(1 + d m>0 ) 

(5) 



where the YJ" l (17) are the spherical harmonics depending 
on the angular part fl of r. In Eq. ([5]) the symmetry of 
the Hamiltonian under reflection at the collision plane 
given by b and v is used in order to reduce the number 
of states. The radial function <Mi( r ) is further expanded 
in terms of (k — l)th order B-spline functions. Converged 
results were found confining the entire space of the elec- 
tron to a sphere of radius r max = 200 with fixed boundary 
conditions. The range [0, r lrl! , x l is divided into N r — 1 in- 



tervals between the knot points (ri = 0, 



, r Nr = 



k+N r 



E 



(6) 



The expansion coefficients a? nl are determined by diago- 
nalizing the atomic Hamiltonian Ho. 

Substitution of the total wavefunction in Eq. ((T|) 
by its expansion given in Eq. ([5]) results in a system of 
coupled equations, 



l dt Cn>V 



'(<) 



E 

nlm 



Cnlrn\J') \ fyu'l'm' I Vi„t I <i> nlm > (7) 

x exp[i(e n 'j/ - e nl )t] . 



These coupled equations are solved for N b fixed val- 
ues of the impact parameter 6 with the initial condition 
Crdm(t = -oo,6) = 5 nlm,mUmi where the index mkmi 
represents the initial state of the atom. The transition 
probability P n im(b) into the atomic state (n, l,m) after 
the collision is given by 



Pnlm(b) = | C n l m {t = +00,6) | 2 



(8) 



The cross section a n im for the transition into the state 
{nlm) follows from 



Onlra — 2 7T J &bb P n l m (b) 

The total cross sections for ionization 
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<=n!m>0 

and for excitation of the target atom 

E ° n 



emit <£ n i m <0 



(9) 



(10) 



(11) 



can readily be calculated where t ni i i is the energy of the 
initial state of the target atom. 
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Figure 1: (Color online) p - Li(2s) collision: The convergence 
behavior at different impact energies is shown for the three 
basis sets A4 (-), A6 (- -) and A8 (• — ■) with maximum 
angular momenta i max = 4, 6 and 8, respectively, (a) To- 
tal excitation probability P(E) EX weighted with the impact 
parameter 6 as a function of b. (b) As (a), but ionization 
probability P(J5) ION . 



energies is investigated. Thereafter, the results for pro- 
ton and antiproton collisions with the alkali metal target 
atoms Li, Na and K are presented and the findings are 
compared with data from literature. Additionally, the 
antiproton cross sections are compared with calculations 
including a hydrogen atom as target. 



Convergence behavior and b-dependent 
transition probabilities 



III. RESULTS AND DISCUSSION 

In what follows, first the convergence behavior and sec- 
ond the dependence of the ionization and excitation cross 
sections on the impact parameter 6 for different impact 



In order to discuss the convergence of the results the 
behavior of the product P(b, E) b is investigated. This 
quantity results after integration over 6 - cf. Eq. ([9]) - 
in the final cross section. Thereby, the probability P 
for a certain transition depends on the impact param- 
eter 6 and the impact energy E. In figures []Ji and [TJs 
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Figure 2: (Color online) p - Li(2s) collision: The convergence 
behavior at different impact energies is shown for the five 
different basis sets P6 (-), P8 (- -), PlOa (• - •), P14 (■ ■ ■ ) and 
PlOb (-) with maximum angular momenta Z max = 6, 8, 10, 14 
and 10, respectively, (a) Total excitation probability P(iJ) EX 
weighted with the impact parameter 6 as a function of b. (b) 
As in (a), but ionization probability P(E) ION . 



the transition probabilities for total excitation and ion- 
ization, respectively, of p - Li(2s) collisions calculated for 
three different energies are presented. The results for p 
- Li(2s) collisions for the same parameters are shown in 
figure [H 

In the case of antiproton collisions with lithium cal- 
culations for three different basis sets A4, A6 and A8 
with maximum angular momenta ^ max = 4, 6 and 8, re- 
spectively, are shown in figure [T] In table [TT] the max- 
imum angular and magnetic quantum numbers as well 
as the total number of basis states are given for the 
basis sets. The calculations using A4 and A6 are fully 
converged. Therefore, the range of integration has been 
changed from -z min = z max = 60 to -z min = z max = 180 
for the calculation using the basis set A8a. This results 
in slightly higher ionization and lower excitation prob- 
abilities at low energies. It was found that the results 



for the calculation using Z max = 8 converge quickly with 
increasing maximum projection of the angular momen- 
tum m max . Consequently, for the subsequent calculations 
dealing with antiprotons as projectile the basis set A8b 
with ^ max = 8, m max = 3 and z max = 180 was chosen 
resulting in a set of 1620 basis functions. 

In the case of proton collisions with lithium calcu- 
lations for five different basis sets P6, P8, PlOa, P14 
and PI 0b different basis sets with maximum angular 
momenta ^ nax = 6, 8, 10, 14 and 10, respectively, are 
shown in figure [2] The further parameters of these ba- 
sis sets are given again in table [ill In contrast to the 
antiproton calculations, much higher angular momenta 
are needed to achieve convergence - especially for low 
energies E < 4keV. The results for the basis sets PlOa 
and P14 both with m max = 3 seem to be converged with 
respect to I. Hence, i max = 10 was chosen but m was in- 
creased to m max = 6 leading to the basis set PlOb. Addi- 
tionally, the integration range z max was also enlarged to 
Zmax = 180. For all considered energies with E > 4keV 
these parameters lead to satisfyingly converged results 
and were therefore used for all proton collision calcula- 
tions. 

From the previous analysis it can be concluded that 
convergence is achieved much faster (i) for antiprotons 
than for protons, faster (ii) for excitation than for ion- 
ization and (iii) at higher than at lower energies. Figures 
[1] and [2] also provide insight into the physics of the col- 
lision process. For high energies the same behavior for 
antiproton and proton collisions can be observed. For 
energies below the validity regime of the first Born ap- 
proximation the transitions in antiproton collisions take 
place at smaller impact parameters compared to protons. 
For close encounters which are more important for low 
energies the advent of the projectile inside the orbit of 
the target electrons creates in the case of protons an in- 
creased or for antiprotons a decreased binding of the elec- 
trons. This situation leads to a decrease (p) or increase 
(p) of P for small b [24| and a shift of the proton P curves 
to larger b. The ionization probability is for p and p more 
concentrated in the vicinity of the nucleus. This can be 
explained using the simple picture that the mean veloc- 
ity of the electrons close to the nucleus is higher than 



basis /max 


W max 


states 


basis 


/max 


771 max 


states 


A4 4 


4 


810 


P6 


6 


6 


1188 


A6 6 


6 


1512 


P8 


8 


8 


1620 


A8a 8 


8 


2430 


PlOa 


10 


3 


2052 


A8b 8 


3 


1620 


P14 


14 


3 


2916 








PlOb 


10 


6 


3024 



Table II: Parameters of the basis sets used for the conver- 
gence studies. Basis sets beginning with A (P) are used in 
calculation with antiprotons (protons). For each basis set the 
maximum angular quantum number i m ax, the maximum mag- 
netic quantum number m max and the total number of basis 
states are given. 



at larger distances and therefore less energy transfer is 
required. On the other hand the excitation probability 
has a longer tail for large b compared to ionization. Par- 
ticularly at high energies care has to be taken that the 
calculations converge in the considered impact parameter 
range. 



B. Cross sections for proton collisions 

Especially for Li atoms but also for Na a number of the- 
oretical and experimental results as well as derived fits 
exist in literature. Thereby these two collision systems 
become good candidates to test the used method. Ad- 
ditionally, the results of different theoretical approaches 
and by that their applicability can be compared. The 
achieved understanding of the proton systems may be 
used for the discussion of the antiproton collisions later 
on for which a thorough comparison is not possible due 
to the sparseness of literature dealing with antiproton 
alkali-metal atom collisions. The present results for pro- 
ton collisions with K complement the sparse literature 
data for this collision system. 



1. Ionization 

In figures [3j |4] and \5\ the results of the calculations for 
proton collisions with Li(2s), Na(3s) and K(4s), respec- 
tively, are presented. The cross sections for the ionization 
of alkali metal atoms A, 



p + A(n l s) 



A 4 



(12) 



where A stands for either Li, Na or K initially in their 
ground states n^s, e.g., Li(2s) are shown in the subfigures 
[3^,, [4^ and [5^, respectively. The ionization cross section 
for proton collisions includes two processes. First the ion- 
ization of the alkali metal atom due to electron capture by 
the proton and second the ionization of the electron into 
the continuum. The sum of both cross sections is some- 
times also referred to as electron loss cross section. The 
electron capture by the projectile is the dominant process 
for low energies but vanishes fast with increasing energies 
and becomes negligible for E > 100 keV. For intermedi- 
ate and high energies the ionization into the continuum is 
the dominant electron loss process. Therefore, in the fol- 
lowing discussion the present ionization calculations are 
also compared with electron capture cross sections from 
the literature at low energies and with literature results 
for ionization excluding electron capture by the proton 
at high energies. 

The present proton ionization cross section for Li(2s) 
in figure [3^ matches perfectly with the 2s contributions 
of the theoretical ionization cross section by McCartney 
et al. Lyi and also with the experimental results by Shah 
et al. [3|. The ionization cross section by Schweitzer et 
al. [l[ is somewhat smaller at high energies. The con- 
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Figure 3: (Color online) p - Li(2s): (a) Ionization and cap- 
ture. Theory (ionization): solid curve, present results; doubly- 
dashed dotted curve, Schweinzer et al. [lj; dashed curve, 
McCartney et al. [18(1 (Li (2s)). Fit(ionization): thin solid 
curve, Wutte et al. 0. Experiment(ionization): squares, 
Shah et al. 0|; crosses, Shah et al. 0] (Li(2s)). Fit(capture): 
dashed doubly-dotted curve, Tabata et al. [J]; dashed dotted 
curve, Morgan et al. ffj. Experiment(capture): diamonds, 
Aumayr et al. jjj. (b) Total excitation and excitation into 
Li(2p). Theory (total excitation): solid curve, present re- 
sults. Theory(excitation into Li(2p)): dotted curve, present 
results; dashed dotted curve, Brandenburger et al. Q; dashed 
doubly-dotted curve, Stary et al. :17j; long-dashed curve, 
MC2 Nagy et al. @|; short-dashed curve, MC3 Nagy et al. 
@|. Fit(Li(2p)): thin solid curve, Wutte et al. @]. Experi- 
ment(Li(2p)): diamonds, Aumayr et al. 



tribution of the Is electrons to the ionization cross sec- 
tion which is not included in the present calculations has 
been determined theoretically by Sahoo et al. [25[ and 
McCartney et al. [l8| and experimentally by Shah et al. 
[3j. For energies smaller than 100 keV the contribution 
of the inner shell becomes negligible compared to the 
one of the outer shell. For high energies the Is contri- 
bution is in accordance with the difference between the 
present 2s results and the Li electron ionization cross 
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Figure 4: (Color online) p - Na(3s): (a) Ionization and cap- 
ture. Theory (ionization): solid curve, present results; long- 
dashed dotted curve, Stary et al. [171 ]: short-dashed dot- 
ted curve, Zapukhlyak et al. Experiment (ionization): 
squares, Zapukhlyak et al. [io|. Fit(capture): dashed curve, 
Tabata et al. 0); doubly-dashed dotted curve, Morgan et al. 
pj. Experiment (capture): diamonds, Aumayr et al. [ill ], (b) 
Total excitation and excitation into Na(3p). Theory (total 
excitation): solid curve, present results; long-dashed curve, 
Shingal et al. p2]. Theory (excitation into Na(3p)): dotted 
curve, present results; dashed doubly-dotted curve, Jain et al. 
[l3| : dashed dotted curve, Shingal et al. doubly-dashed 
dotted curve, Shingal et al. [l2| ]; short-dashed curve, Stary et 
al. [T3|. Experiment(Na(3p)): diamonds, Aumayr et al. [llj ]. 



section. For energies smaller than lOkeV the electron 
capture by the proton becomes the dominant ionization 
process. Down to 4keV the present findings are in good 
agreement with literature results shown in figure [3^ for 
capture and ionization. However, for energies smaller 
than 4 keV the present ionization cross section is clearly 
smaller than all other shown results. This is in accor- 
dance with the difficulty to achieve convergence in the 
energy range E < 4keV for proton collision already dis- 
cussed in section IIII Al This may be a consequence of 
the one-center approach which is expected to converge 



Figure 5: (Color online) p - K(4s): (a) Ionization and 
capture. Theory (ionization): solid curve, present results. 
Theory (capture): long-dashed dotted curve, Fritsch fl5l ] . 
Fit (capture): dashed curve, Tabata et al. [4J; dashed dot- 
ted curve, Morgan et al. [f|. Experiment(capture): squares, 
Gieler et al. (b) Total excitation and excitation into 

K(4p). Theory(total excitation): solid curve, present results. 
Theory(excitation into K(4p)): dotted curve, present results. 
Experiment (K(4p)): squares, Gieler et al. [la ]. 



slowly when trying to describe the electron capture pro- 
cess properly. 

In figure 0^ the results for the ionization of the Na 
atom initially in the ground state are shown. The find- 
ings are in good agreement with the recent results by 
Zapukhlyak et al. [101 ] . Here, especially their experimen- 
tally values match with the present curve except for the 
three last data points with E > 17keV which show an 
unexpected behavior. Their theoretical ionization cross 
section agrees for energies higher than 6 keV with the 
present one but is larger for smaller energies. The cross 
section of Stary et al. [I?) which also covers the range 
from low to high energies differs from the present find- 
ings as well as from the literature results shown here. 
The data for electron capture by the proton shown here 
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are consistent with the present findings for the ioniza- 
tion cross section. However, the maximum of the fit by 
Tabata et al. Q has a clearly higher value. Again, it 
can be observed that the present ionization cross section 
for proton collisions is not fully converged for E < 4 keV 
leading to smaller values in this energy range. 

In figure^ the results of the proton - K(4s) collision 
calculations are presented. For potassium targets the lit- 
erature data on proton cross sections are sparse. To the 
best of the authors' knowledge no ionization cross sec- 
tions for proton collisions exists. Therefore, the present 
ionization cross section may be compared with results for 
electron capture. However, this is only meaningful for 
low energies E < lOkcV where electron capture is the 
dominant ionization process. The calculations by Fritsch 
[lH , the experimental data measured by Gieler et al. [l6[ 
as well as the fit provided by Tabata et al. [4[ of the 
electron capture cross section are in accordance with the 
present ionization cross section for E > 4keV. The fit- 
ted capture cross section by Morgan et al. [|[ results in 
lower values in the relevant energy range between 4 and 
lOkeV. In the high-energy regime the present cross sec- 
tion shows the same qualitative behavior which already 
has been observed for Li and Na. 



2. Excitation 

In the subfigures [3)3, 0J) and the proton excitation 
cross sections for Li(2s), Na(3s) and K(4s) are shown. 
The total excitation of an alkali metal atom A initially 
in its ground state ri^s, 

p + A(n;s) -> p + A(nZ), (13) 

which is the sum of all cross sections for transitions into 
excited states nl, cf. Eq. (fTTj) , is given. Additionally, 
the cross section for the excitation process into the first 
excited state n^p of A, 

p + A(n !: s) -» p + A(n iP ), (14) 

is given, too. The excitation into the first excited state 
riip is the dominant excitation process especially at high 
energies. Therefore, there are experimental data for this 
excitation transition. It was found in the present inves- 
tigation that it is essential in particular for high energies 
to extend the range of the impact parameter b to values 
up to 90 a.u. in order to achieve excitation cross sec- 
tions which are converged with respect to b. The curves 
in figures Q] and [5] for 500 keV already indicate that the 
transition probabilities for excitation vanish slowly with 
increasing b. 

The present proton excitation cross sections for Li(2s) 
are shown in figure[3jD. To the best of the authors' knowl- 
edge for p - Li collisions there arc no data to compare the 
present total excitation cross section with. For the excita- 
tion into Li(2p) the present results are in good agreement 
with literature data also shown in figure [3Jd, except with 



the calculations by Stary et al. [I?) • Their findings differ 
for E < 6keV and E > 100 keV from all results shown 
here. The experimental data by Aumayr [§] lie for all en- 
ergies below the present calculations. On the other hand, 
the calculations by Brandenburger et al. [7[ and also the 
fit provided by Wutte et al. @ match with the present 
data in the whole energy range. Wutte et al. based their 
fit in the high-energy range on scaled-electron-impact ex- 
citation cross sections. The calculations by Nagy et al. § 
were performed with multi-configuration wavefunctions 
with an orbital basis up to n — 2 (MC2) and up to n = 3 
(MC3). 

In figure [4}d p - Na(3s) collision cross sections for the 
total excitation and excitation into the 3p state of the 
sodium atom initially in the ground state are shown. 
The theoretical data for the total excitation cross sec- 
tion by Shingal et al. [T4] agree well with the present 
results although they show a feature around 4 keV which 
is not reproduced by the present findings. However, their 
excitation cross section into the Na(3p) state follows al- 
most completely the present results - even around 4 keV. 
The older calculation by Shingal et al. [3] agrees rea- 
sonable in the energy range 4 — 14 keV but shows a dif- 
ferent behavior for higher and lower energies. Although 
the calculations of Jain et al. [l3| lead for all energies to 
higher values their qualitative behavior is comparable to 
the present results. The findings of Stary et al. [ItJ show 
the same behavior as their results for p - Li collisions in 
figure [3Jd, namely, a cross section which is comparable 
around the maximum but falls off too fast for higher and 
lower energies. The experimental data provided by Au- 
mayr et al. [Til ] is in line with the present cross section 
for excitation into Na(3p). It also shows a feature around 
4keV. 

In figure [5] the results of the present p - K(4s) collision 
calculations are presented. For excitation into the K(4p) 
state the experimental findings of Gieler et al. [IB] are 
in good agreement with the present calculations around 
the maximum but then start to differ for E < 4keV. 
Their data points fall off faster while the present result 
shows a behavior which has been already observed for p 

- Na collisions in figure UJd. For Na the slope of the curve 
changes characteristically around E = 4keV. However, 
there is no comparable feature for p - Li(2s) collisions. 
Although the excitation results - in contrast to ionization 

- for Li and Na collision seem to be reasonable also for 
low energies it is not possible to quantify how reliable 
the p - K(4s) excitation cross sections for E < 4keV 
are. The splitting of the energy levels due to spin-orbit 
coupling which is neglected in the present investigation 
is supposed to be most relevant for the 4p state of K. 
However, the good agreement of the present results with 
the experimental data by Gieler et al. for E > 4keV 
suggests that the effect due to spin-orbit coupling does 
not play a major role with respect to the level of accuracy 
which is achieved by the present method. 

In conclusion, the comparison of the present proton 
ionization and excitation cross sections with literature 
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Figure 6: (Color online) p - Li(2s): (a) Ionization. The- 
ory: solid curve, present results; dashed dotted curve, Stary 
et al. [T3; dashed curve, McCartney et al. Experi- 
ment: crosses (p-Li(2s)), Shah et al. @|; squares (p-Li(2s) 
and p-Li(ls)), Shah et al. @|. (b) Total excitation and ex- 
citation into Li(2p). Theory (total excitation): solid curve, 
present results. Theory(excitation into Li(2p)): dotted curve, 
present results; dashed dotted curve, Stary et al. [13]. Exper- 
iment (excitation into Li(2p)): diamonds (p-Li(2s)). Aumayr 
et al. %. 



data results in a good overall agreement in the energy 
range 4keV < E < 1000 keV. Thereby, the applicability 
of the present method is confirmed. The findings by Stary 
et al. [17|, however, differ from the present and the shown 
literature results. The calculations for p - K collisions 
complement the data provided by the sparse literature 
on this collision system. 



C. Cross sections for antiproton collisions 

Only very few data for antiproton-alkali metal atom 
collisions exist in the literature. Cross sections are avail- 
able for the ionization of Li(2s) and Na(3s) as well as for 



Figure 7: (Color online) p - Na(3s) (a) Ionization. Theory: 
solid curve, present results; dashed curve, Stary et al. [171 ]. 
(b) Total excitation and excitation into Na(3p). Theory (total 
excitation): solid curve, present results. Theory (excitation 
into Na(3p)): dotted curve, present results. 



excitation into Li(2p) and Na(3p) by Stary et al. [T7| . 
Furthermore, there are ionization cross sections for p - 
Li collisions calculated by McCartney et al. [l|[. How- 
ever, no cross section exists for ionization or excitation 
into K(4p) for K targets and also for all three considered 
target atoms there are no total excitation cross sections 
to which the present results could be compared. 



1. Ionization and excitation 

The ionization cross sections for antiproton collisions 
with the target atoms Li(2s), Na(3s) and K(4s) arc shown 
in figures [S^, UJa, and [5Ji, respectively, and for excitation 
accordingly in figures [Hp, [7b and [5b. The theoretical re- 
sults for ionization in p - Li collisions by McCartney et 
al. [181 ] agree well with the present findings. However, 
they only cover the high-energy regime E > 30 keV. The 
calculated antiproton ionization cross sections for Li and 
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Figure 8: (Color online) p - K(4s) (a) Ionization. Theory: 
solid curve, present results, (b) Total excitation and excita- 
tion into K(4p). Theory(total excitation) : solid curve, present 
results. Theory(excitation into K(4p)): dotted curve, present 
results. 



Na targets by Stary et al. [l7| both differ considerably 
from the present findings. Their p - Li ionization cross 
section behaves differently for low to intermediate ener- 
gies but seems to converge to the present findings for high 
energies. Their p - Na ionization cross section, however, 
shows a different behavior compared to the present curve 
in the whole energy range. 

The cross sections for the excitation into the first ex- 
cited state for Li and Na target atoms calculated by Stary 
et al. *17J both share the same features. Their cross sec- 
tions agree with the present curves around the maxima at 
E w lOkeV and E ss 15keV for Li and Na, respectively, 
but fall off faster for lower and higher energies. The same 
behavior has been observed in the case of proton col- 
lisions in figures [3Jd and HJ) for excitations into Li(2p) 
and Na(3p), respectively. Therefore, their results differ 
once more from the outcome of the present investigation. 
The aim of Stary et al. was to obtain results compara- 
ble to literature data but using smaller basis sets within 



an optical potential approach adapted to this problem. A 
Feshbach projector formalism for the solution of the time- 
dependent Schrodinger equation leading to a finite set of 
coupled-channel equations with complex potentials was 
used. Thereby, two conditions were assumed to be ful- 
filled. First, the interactions occur instantaneously and 
second, the energy distribution of the Q-space which is 
the complement of the finite model space has a peak lead- 
ing to the assumption of an average Q-space energy e. 
Furthermore, a scaling factor is used which restores the 
correct energy dependence of the optical potential and 
which is determined at high impact energies. Since the 
present results for the proton case seem to be more reli- 
able than their one-center calculations the present results 
for antiproton collisions with Li and Na are also consid- 
ered to be more reliable than theirs. If their solutions are 
converged as it was claimed by Stary et al. then either 
not both of the above mentioned conditions are fulfilled 
or the introduced scaling factor has a different functional 
behavior. 

To the best of the authors' knowledge no literature on 
p - K cross sections exist for the considered energy range. 
The present cross sections for excitation and ionization 
of K in figure [5] show a qualitatively similar behavior like 
for p - Na collisions in figure [7] but with higher values 
throughout the energy range. 

Until now, experimental results for the antiproton- 
alkali metal atom collision systems are completely miss- 
ing in the considered energy range. It is remarkable that 
the experimental data of Aumayr Q for excitation into 
Li(2p) by proton collisions fits better to the present an- 
tiproton than proton data. 



2. Comparison of antiproton with proton cross sections 

While for sufficiently high energies a similar behavior 
for proton and antiproton cross sections is to be expected 
the collision processes should differ for lower energies. In 
contrary to the proton collisions no electron capture by 
the projectile is possible for antiprotons. Since the elec- 
tron capture is the dominant ionization channel for low 
energy proton collisions noticeable differences especially 
for the antiproton ionization cross sections can be ex- 
pected in the low-energy regime. In what follows the an- 
tiproton and proton cross sections are compared in some 
detail for high, intermediate and low impact energies. In 
figure [H the ratios of proton to antiproton cross sections 
for ionization and excitation is given for the three con- 
sidered target atoms. 

To begin with the comparison focuses on the high en- 
ergy behavior of the antiproton and proton cross sec- 
tions. In the validity range of the first Born approx- 
imation no differences in the cross sections for differ- 
ent projectiles like electrons, protons and antiprotons 
with the same velocity are expected because in this ap- 
proximation the cross sections only depend on the ab- 
solute value of the projectile charge. It is a high en- 
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ergy approximation. A linear decrease of the ionization 
cross section for all three alkali metal atoms can be ob- 
served on a doubly-logarithmic scale for high energies 
100 keV < E < 1000 keV for protons as well as for an- 
tiprotons. Therefore, the general fit formula 



(£■) = < T i,,„(£V) ( 



(15) 



for ionization cross sections in this energy range can 
be proposed, where (Ji on (Eo) is the ionization cross sec- 
tion for an arbitrary Eq in the range 100 keV < Eq < 
1000 keV and a is a fit parameter which gives the slope 
of the linear curve on a doubly-logarithmic scale. The fit 
parameters which may be proposed for the three alkali 
metals colliding with protons are given in table IIIII The 
fits for Na and K reveal a direct proportionality between 
the ionization cross section and the inverse of the energy, 



1 (E ) Eq 

E 



oc 



1 

E' 



(16) 



in the considered high energy regime. This proportional- 
ity holds approximately also for the present Li ionization 
cross section. The proposed fits which are also shown 
in figures [3^, HJi and [5^ match well with the calculated 
ionization cross sections for E > 150 keV. These fits ob- 
tained for the proton case are also shown in figures [6^i, 
[7£i and [5^ m order to compare them with the antiproton 
ionization cross sections. It can be seen that the proton 
fits match remarkably well with the antiproton ioniza- 
tion results for energies higher than 150 keV. Therefore, 
the antiproton ionization cross sections also decrease pro- 
portional to E^ 1 where this proportionality again holds 
only approximately for Li targets. This means that for 
energies higher than 150 keV no specific features are ex- 
pected for antiproton ionization cross sections with the 
considered alkali metal atoms. And in turn for these ener- 
gies the treatment of proton collisions should be sufficient 
which is especially in the case of experimental studies less 
demanding. 

However, both systems differ strongly regarding ion- 
ization for energies lower than 100 keV due to the elec- 
tron capture process which is only possible for protons. 
The proton cross section is strongly enhanced which can 
also be seen in figure The maxima of the proton and 
antiproton ionization cross sections approximately at 45 



Atom 


Eo 


O-ion(-Eo) 


a 


Li(2s) 


141.3 


2 


-0.9386 


Na(3s) 


138 


2 


-1 


K (4s) 


151 


2 


-1 



Table III: Parameters for the description of the ionization 
cross section for the energy range 100 keV < -Eo < 1000 keV 
using the fit formula (|15[1 . a is a dimensionless fit parameter. 
The energy Eo is given in keV and the cross section o"i on (-Bo) 
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Figure 9: (Color online) Ratio of proton to antiproton cross 
sections <j p /<j ap . The ratios of the ionization and excitation 
cross sections for the target atoms Li(2s), Na(3s) and K(4s) 
are given. 



and 10 [in 10 _16 cm 2 ], respectively, for Li targets differ 
by a factor 4.5. The ionization maxima for Na and K 
targets differ approximately by a factor 5.5 and 6.5, re- 
spectively. The maxima are all located between 4 and 
6keV where the proton maxima tend to occur at lower 
energies than the corresponding antiproton maxima. 

A comparison of the present excitation cross sections 
for proton and antiproton collisions yields that they also 
agree for high energies E > 150 keV. The antiproton 
maximum for Li targets lies around lOkeV and is 10% 
lower than for proton collisions. The antiproton maxima 
for Na and K are situated at approximately 15 keV with 
« 20% smaller values than for the proton case. But below 
their maxima the Na and K excitation curves for antipro- 
tons and protons excitation cross section have compara- 
ble values. 

The most striking feature of figure [5] is that the ratios 
of the proton to antiproton ionization cross sections in- 
crease strongly for low-energy collisions while the ratios 
for excitation only vary comparably weakly around 1. In 
the case of ionization the electron capture channel be- 
comes important for low-energy proton collisions leading 
to large ionization cross sections compared to antiproton 
collisions. In the case of excitation for both projectiles 
the same channels are open. 

It can be concluded for the antiproton cross sections 
that the present results complement and improve the ex- 
isting data on antiproton-alkali metal atom collisions. 
While the excitation cross sections are comparable for 
proton and antiproton projectiles the proton ionization 
cross sections are strongly enhanced at low energies due 
to electron capture. For high energies E > 150 keV pro- 
ton and antiproton collisions with Li, Na and K result 
in the same ionization cross sections which decrease pro- 
portional to E~ l . 
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values compared to the proton collision systems. 

In figure [TUb the total excitation cross sections for the 
three alkali metals are compared. Although the overall 
behavior of their excitation cross sections are similar they 
differ in detail. On an absolute scale the values of the 
antiproton excitation maxima for the three atoms differ 
considerably [in 10" 16 cm 2 ]: Li = 46.2, Na = 39.7, K = 
65.9, i.e., the maximum for K is 66% higher than that for 
Na. This can be made plausible by comparing the energy 
differences of the ground states to the first excited states: 
Li = 0.068, Na = 0.089, K = 0.059 since the n,s -> n lP 
transition is the dominant excitation channel. 

The cross sections for ionization and excitation of the 
hydrogen atom by antiproton impact were also calculated 
which are presented in figure 1101 too. The qualitative 
behavior of the hydrogen cross sections is comparable 
to those of the alkali metal atoms reflecting the shell 
structure of the alkali metal atoms with an outer va- 
lence electron in an s state. However, the absolute values 
of the cross sections differ clearly. First, the cross sec- 
tions for hydrogen are much smaller due to the tighter 
binding of the electron which leads to a higher ioniza- 
tion energy and a smaller spacial extension. Second, the 
maxima are shifted to higher impact energies which can 
be explained by the higher average velocity of the elec- 
tron in the ground state of hydrogen. For high energies 
(E > 1000 keV) the ionization cross section of hydrogen 
seem to approach those of the alkali metal atoms. At 
these energies the ionization cross section of hydrogen is 
expected to decrease like E^ 1 \ogE. 



Figure 10: (Color online) (a) Single-electron ionization cross 
section for antiprotons colliding with Li(2s), solid curve, 
Na(3s), dashed dotted curve, K(4s), dashed curve, and H(ls), 
dashed doubly-dotted curve. Additionally, the fits describ- 
ing the high-energy behavior of proton scattering are given, 
(b) Total excitation cross section for antiprotons colliding 
with Li(2s), solid curve, Na(3s), dashed dotted curve, K(4s), 
dashed curve, and H(ls), dashed doubly-dotted curve. 



D. Comparison of the antiproton cross sections 

In figure [TUb the ionization cross sections for the three 
alkali metal atoms Li, Na and K colliding with antipro- 
tons are plotted together with the high-energy fits ex- 
tracted earlier from the proton calculations. The qualita- 
tive behavior of the cross sections for these atoms is simi- 
lar in the whole energy range. All curves converge to their 
corresponding proton results for energies E > 150 keV 
which lie close to each other. The differences between the 
heights of the ionization cross sections for E < 100 keV 
as well as the ordering of the curves may be explained by 
the different ionization energies: Li = 0.198, Na = 0.189, 
K = 0.160 which can also be found in table|TJ All maxima 
lie around 4 to 6 keV which is somewhat below the aver- 
age velocity of the valence electrons and have far lower 



IV. CONCLUSION 

Time-dependent close-coupling calculations of ioniza- 
tion and excitation cross sections for antiproton and pro- 
ton collisions with the alkali metal atoms Li(2s), Na(3s) 
and K(4s) have been performed in a wide energy range 
from 0.25 to 1000 keV. The target atoms are treated as 
effective one-electron atoms using a model potential. The 
total wave function is expanded in an one-center ap- 
proach in eigenfunctions of the one-electron model Hamil- 
tonian of the target atom. The radial part of the basis 
functions is expanded in B-spline functions and the an- 
gular part in a symmetry-adapted sum of spherical har- 
monics. The collision process is described in the classi- 
cal trajectory approximation. In the present calculations 
the results converged faster for collisions involving an- 
tiprotons than protons, faster for high than for low ener- 
gies and faster for excitation than for ionization. Good 
agreement with literature data has been achieved for the 
proton-alkali metal atom cross sections for E > 4keV. 
However, for antiproton-alkali metal atom collisions lit- 
erature data are sparse. The comparison to the calcula- 
tions of antiproton collisions with Li and Na by Stary et 
al. shows the same disagreement with the present findings 
as it was found for their proton collision results. In view 
of this disagreement with literature data for proton colli- 
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sions it can be stated that the calculations by Stary et al. 
were cither not fully converged or the assumed conditions 
not fulfilled. To the best of the authors' knowledge the 
first cross sections for p - K collisions in the considered 
energy range are presented. The ionization cross sec- 
tions for protons and antiprotons differ considerably for 
energies smaller than 100 keV due to the electron capture 
process which is only possible for protons and is the dom- 
inant ionization channel at low energies. The qualitative 
behavior of the antiproton cross sections is comparable 
for all three alkali metal atoms but differs in the abso- 
lute values depending on the atom-specific ionization and 
excitation energies. A comparison with hydrogen as tar- 
get atom yields the same characteristics as for the alkali 
metals due to the common s state structure. However, 



the cross sections of hydrogen have much lower values 
and the hydrogen ionization and excitation maxima are 
shifted to higher impact energies because of the tightly 
bound Is electron. For the proton ionization cross sec- 
tions a simple fit formula is proposed for the energy range 
from 150 to 1000 keV which also describes the properties 
of the antiproton ionization cross sections in this energy 
range well. The fit reveals that the ionization cross sec- 
tions decrease proportional to E^ 1 in this energy range. 
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